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<^ : Abstract 

Tij- I Using a recent calculation of the perturbative hard part for dilepton production that sums 

^ I large threshold corrections to all orders in perturbative QCD, we compute the corresponding 

cross sections. The hard part has been evaluated using principal value resummation and 
contains all singular momentum-dependent corrections. We also include a resummation of 

Q^ . large Sudakov terms, which are independent of parton momenta. We give predictions for the 



dilepton-mass distribution, the rapidity distribution and the rapidity-integrated K-factor 

CL|| at fixed-target energies and compare with various experimental results in several kinematic 

p • regimes. We find that principal value resummation produces cross sections that are finite 

^ ■ and well-behaved. For both protons and anti-protons on fixed targets, the resummed cross 

sections are, in general, in excellent agreement with the data. 
^ ■ 

X 
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1 Introduction 

It is well-known that certain liadron-liadron scattering cross sections, evaluated in perturba- 
tive QCD (pQCD), suffer from the presence of large corrections at all but very high energies 
[l[|. These are associated with soft gluons that produce large threshold effects. The primary 
example of such processes is the dilepton-production cross section, in which the large cor- 
rections are associated with the region where all the available partonic momentum is carried 
away by the dilepton final state. We recently completed the calculation of all these large 
corrections to all orders in pQCD based on the method of principal value resummation 
0, whose main advantage is that it by-passes the infrared (IR) singularities associated with 
the Landau pole of the QCD running coupling Q . The corresponding resummed hard part 
is thus independent of arbitrary IR cutoffs, and can be used to produce unambiguous pre- 
dictions for the cross section. These are, moreover, stable from a perturbative point of view, 
since they include all such large corrections. 

The principal value prescription exponentiates the corrections into a leading and next-to- 
leading exponent, E^ and E^l, respectively. These exponents contain the effects of the one- 
{El) and two-loop {E^l) running coupling and asymptotically reproduce the resummed 
perturbative series in the partonic phase space region where perturbation theory is valid. 
Therefore, the corresponding resummed hard part is well-defined mathematically through- 
out the whole kinematic region, with all the correct properties associated with perturbation 
theory. Furthermore, to completely calculate these exponents, which resum the large thresh- 
old corrections to all orders in perturbation theory, only one- and two-loop calculations of 
the corresponding hard part are needed. In [Q and 0, we have performed an exhaustive 
study of the properties of these exponents and the associated resummed hard part in the 
context of the dilepton cross section. The main objective of the present work is to use this 
hard part to compare predictions for the physical cross section and related quantities with 
experiment. 

In section 2 we will review the general formulation of the dilepton production cross sec- 
tion, both with the resummed hard part and with finite-order perturbative expressions up 
to two loops p . Our formalism is valid for the differential cross section, with respect to the 
dilepton invariant mass, as well as for its derivative with respect to rapidity, y, a.t y = 0. In 
section 3 we review and discuss the resummation of large Sudakov constants , which are 
particularly important as well at low invariant mass. In section 4 we explicitly exhibit the 
predictions resulting from our formalism, i.e., mass and mass-rapidity (or mass- a; i?) distribu- 
tions and i^-factors, and compare these with data from available fixed-target experiments. 
We will find that principal value resummation produces corrections that are finite and not 
very large compared to 2-loop corrections. As one goes to higher mass values, the gap be- 
tween the resummed and finite-order cross sections diminishes. We will also see that the 
resummation is sensitive to the choice of parton distributions, with global distributions lead- 
ing to larger cross sections than those obtained from older deep inelastic scattering (DIS) 
fits. Finally, in section 5 we summarize and discuss our findings. Some numerical details 
are discussed in an appendix. We postpone comparison with collider experiments to future 
work. 



2 The inclusive cross section for dilepton production 

In this section, we review the general formulation for the cross section for dilepton produc- 
tion in hadronic collisions and provide all the necessary formulas for our calculations. This 
will streamline the discussion which we will give below and facilitate comparisons with var- 
ious experimental results. For completeness, we present the general formalism for the cross 
section, valid at either fixed-target or collider energies, even though in this paper we will 
deal only with the former, i.e., for invariant masses much smaller than the W or Z mass. In 
the following, we will follow the notation of 0, |^. In addition, we will also use the notation 
that for any physical quantity of interest R, Rl stands for this quantity calculated using 
principal value resummation with the leading exponent E^ only, R for the same resummed 
quantity using the full exponent E = El + Ej^l, while in finite-order pQCD the same quan- 
tity will be written as Ri = Y^]=o a^R^-'\ where R^^^ is the j-th order perturbative correction 
and a = as{Q'^)/'K. 
In the reaction 

K{pi)+h2{p2)-^lKQn+X, (1) 

the dilepton pair may be produced by a neutral virtual gauge boson V {V E {7, Z}). We 
will use the standard notation s = {pi +^2)^7 t = Q'^/s, where Q is the dilepton invariant 
mass. 

2.1 Mass distributions and iC- factors 

We will first deal with the single-differential inclusive cross section 

da"" 



dQ' 



aUQ')W''{T,Q'), (2) 



where a^ is the point-like Born cross section for qq ^ V* -^ II and contains all the elec- 
troweak parameters and physical dimensions. In this paper, we are mostly concerned with 
large perturbative corrections which come from the diagonal-flavor quark-antiquark sub- 
process. However, when we actually calculate the cross sections, we will also include the 
contributions coming from the non-singular terms in the 1-loop hard part. Furthermore, we 
will use parton distribution functions normalized to deep-inelastic scattering experimental 
data (DIS scheme). 

We may now write the hadronic factor as 

tt ^0 Xi X2 
= Y^ClfdzuM^a)^^^^^^^^, (3) 

. J T Z 

a,b 

where the sum extends over all active parton (quark, antiquark and gluon) flavors, 

n'y — n'y — n^ — n^ — ^2 

C^j = Cf^ = C^j^Cf = l + il-A\eMn'ew)^ (4) 



where e/ is the fractional charge of quark / and 6w is the electroweak mixing angle. Finally, 
the parton luminosity is defined as 

/* J- /7 T* /f'^y 

J^ab{,r/z)= 5{1-t/{zXiX2)) ^ -Fa/hi{Xl,Q'^)Fb/h2{x2,Q'^)- (5) 

Jo Xi X2 

For our purposes, it is convenient to integrate eq. (^ by parts to obtain the equivalent 
expression 

For the purely electromagnetic case, ^ = 7, we have 

where N is the number of colors and a^ is the (electromagnetic) fine structure constant. 
Correspondingly, for V = Z we have 






^ Z^ll 



(Q2 _ MIY + M|r| 



ANsin^ewcos^OwMz 
where Mz = Z boson mass, T^^ii is the partial decay width of Z to one dilepton species, 

aeM^(l + (l-4sin2^vi/)^ 



z^a 



48 sin^ Ow cos^ 6 



(9) 



w 



and Tz is the total width of the Z (summed over all leptonic and hadronic decay channels). 

The physical cross section da/dQ"^ will in general involve the sum of da^ jdQ^ over both 

vector bosons, as well as an interference term, proportional to cr^ , between the amplitudes 

mediated by a photon and a Z 0. For dilepton production in the continuum, i.e., away 

from resonances, which at fixed-target energies i^ Q^ <C M|, the quantities af, a^ are 

suppressed relative to a^ by {Q /MzY and {Q /MzY, respectively, and therefore we may 

write 

da da'^ 

dQ^ ~ dQ^' ^ ' 

The i^-factor is a quantity that nicely measures the size of the radiative corrections in 
the hard part by more-or-less cancelling the effect of non-perturbative parton distributions. 
It may be defined as 

_ da/dQ^ 

dao/dQ^' ^ ^ 

where the denominator is calculated in the same way as the numerator, eqs. (0), (0), but 

using only the tree-level hard part. In the continuum, and in accordance with eq. (p!0D , we 

may write 

daVdQ^^W^ir^ 



^This is the region where aU fixed-target experiments are performed since for these experiments, s < 
M|. For simphcity, we will exclude hadronic resonances, like J/^ and T, which are subtracted from most 
experimental data as well. 



Since we will only be interested in fixed-target energies, these are the formulas we will use 
throughout this paper. 

We may now obtain resummed expressions for the above quantities by using the re- 
summed hard part ujab in the above formulas. The expression for the resummed hard part 
is 0, |, ID u)ab{z, a) = SafS^fUJffiz, a) with 



'ff^ 



ujff{z,a) = A{a)I{z,a), 



(13) 



where the function A{a) exponentiates the (5(1 — z) terms and will be analysed in detail 
in section 3, while I{z, a) exponentiates the momentum-dependent plus-distributions that 
produce large perturbative corrections at the edge of phase space z ^>- 1. The form of 
I{z, a) depends on the finite-order structure of the plus-distributions to be resummed and, 
for best accuracy, should reproduce that structure upon a Taylor expansion in a. In [|| we 
generically studied I{z, a) using a simplified expression that reproduces most of the finite- 
order large threshold corrections up to two loops and resums the bulk of these corrections to 
all orders. That is sufficient for examining all the properties of the resummation. Since in 
this work we will compare with experiment, however, we will use a more accurate expression 
for the hard part, that reproduces all of the finite-order large threshold corrections up to 
two loops and resums them to all orders. Additional resummed structures in I{z,a), that 
are calculable in principal value resummation ^, |^ and would reproduce, upon expansion, 
subdominant remaining threshold corrections starting at three-loop order, are negligible. For 
completeness, we reproduce the expression for I{z, a) below: 



I{z,a) =6{l-z) 



-r(l + Pi)Q[Pi,P2] 



(14) 



+ 



with 



Q[Pl,P2] 



sin(7rPi) 



TT 



l + P2[^\l + Pl)+^^'\l + Pl] 



-TX 



-2P2COs(7rPi)^(l + Pi), (15) 



where E = El + Enl contains all the large logarithmic corrections, F is the Euler Gamma 
function and \1/, '^^^'' are the usual polygamma functions, and 



Pk{x,a) = 



d^E{x,a) 
k\dx^ 



(16) 



where it is understood that E = E{x, a) is a function that has a polynomial representation 
in a; = ln(l/(l — z)) and a in the perturbative regime. From eqs. (^ 



I, (|T3D and (|T|), we 

may eliminate the (5-function and plus-distributions by explicitly integrating the former and 
using the property J^ dz'[f{z')]^ = — Jq dz'f{z') for the latter. Hence, our final form for the 
resummed hadronic function is 



W''{T,Q')=J2CjA{a) / dz[l 
f ■'^ 

where the function 7Y(z, a) is given by 



1-L{z,a)] 



dz \ z 



(17) 



H{z,a) 



rfxe^^^'^^F (1 + Pi (x, a)) Q[Pi{x, a), P2(x, a)] 



;i8) 



It will also be of interest to compare the resummed expression for the cross section or 
the /^-factor with the corresponding quantities obtained from finite-order calculations. The 
corresponding hadronic factor is given by 



d f Tab{T/z) 

dz\ z 



Wr{r,Q^)=Y.Cl fdz\fdz'uU^\a) 

, Jt \_J z 

a,b 

where the finite-order perturbative hard part is denoted by 

i 

The corresponding finite-order mass-distribution and i^-factor are given by 
-— = a^iQ )Wi {t,Q ) and K^ - 



(19) 



(20) 



W^{t,Q^) 



(21) 



The hard part cj* ? has been calculated up to second order in both DIS and MS schemes, 
[0]-[|TT|, and 0, respectively. The full two-loop corrections are quite lengthy but the 
leading corrections, as in the one-loop case |T^, are the "singular" ones in the z ^ 1 limit. 



m 



i.e., the ones that are proportional to a delta function or a plus distribution. These come 
from virtual and soft-gluon corrections and are only present in the qq channel of the initial 
partons. We reproduce these corrections in the DIS scheme below, together with the leading 
order hard part: 



UJ 



(0) 
ab 



SafSbf-Sil - Z) 



UJ 



(1) 



Cp_ 
2 



2(1 + z^)V^{z) + 31)0(2) + (1 + ^7r2)5(l - z) - 4z - e") 



^f9 






Cf3 

~Y8 



z' + (l- zf 
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ln(l -z) + -z^ -^z + - 



(22) 
(23) 

(24) 



UJ 



(2) 



^\cl{d.2V^{z) + 12V2{z) + [64C(2) + ^2\Di{z) + [112C(3) + 24C(2) + l^V^{z] 



338 



9 

+nfCFW2{z) - ^Vi{z 
"548 



16C(2) 



+CACi.(-44p2(^) + 
44 
"9 

C(2)2 + 120C(3)-3C(2) 



10 + yC(2) 



V,{z) + 



57+yC(2)-24C(3) 



^0(2; 



+5{l-z)[Cl 

+CaCp 

+nfCF 



215 2098,, , 196,, , 154,, ,„ 

- + — «2' - -«3) - -«2' 

16,, , 340,, , 38' 
yC(3) - -C(2) - ^ 



where ({s) is the Riemann Zeta function and 



(25) 



(26) 



The ~ in eq. (p5[ ) signifies that we have dropped all terms that are non-singular in the 
z — i> 1-limit. This, of course, raises the question of how dominant the above singular terms 
are relative to the non-singular ones, that occur both in the qq channel and in the other initial- 
parton-state channels, or conversely, how important are the contributions due to the non- 
singular terms. A thorough numerical investigation of this issue has already been conducted 
0, |10|, |Tl[|. The conclusions vary depending somewhat on the factorization scheme, and 
on the range of energy and invariant mass.0 We may say, in general, that at fixed-target 
energies, the 2-loop non-singular terms account for no more than a 10% negative contribution 
to the exact mass distribution over all channels, with the approximation improving as r 
increases. For details regarding variations of this result with the order, the factorization 
scheme, etc., see |jll|. Relative to the resummed cross section then, the 2-loop non-singular 
terms would only lead to a 3 to 7% reduction. Given these general remarks, and because of 
the complicated structure of the 2-loop non-singular contributions, we did not include them 
in our cross section calculations but we should bear in mind the reduction that they would 
induce in our resummed predictions. As a last note, the terms proportional to 5(1 — z) in 
eqs. ( PBD and (pSj), will also be useful in determining the factor Aipi) in eq. (|r^), as will be 
shown in section 3. 



2.2 Differential mass distributions 

We can also use principal value resummation to compute the rapidity distribution of dileptons 
in the central rapidity region. This is possible because the rapidity, or any other smooth 
function of the vector boson momentum, when in the central region and expressed in terms 
of the initial- and final-state momenta, is approximated by a function of the initial parton 



momenta with corrections that are suppressed by vanishing final-state phase space |12[. To be 
specific, it has been shown in |T^ that for any kinematic variable f] which is a smooth function 
of the vector boson momentum, the resummed quantity (Pa^ /dQ'^df] is approximated by 



2^V 



d'a 



dQ'^df] 



r. 



,Q\n) 



<^'fl 



1 dxi dx2 

Xi X2 



ujff{T/xiX2, a)Ff/h^ (xi, Q'^)Fj/j,^ {x2, Q"^ 



X S{r]-r]{xipi,X2P2)) = 



d^gV 
dQ'^dfj 



[r,Q',v) 



(27) 



where u 



ff 



z, a) is given by eq. (|13|) and ri{xipi, X2P2) denotes the corresponding parton-model 

as 



expression. By analogy to eqs. (|]) and @, we can re- write eq. f|2 



where 



dQ'^drj 



rlj{r/z) 



[r,Q\v) 



2 ^^ ^ ^v 



f 



dz 



dz'ujjj'{z', a) 



d_/y^1f{r/z) 
dz\ 



ff^ 



{2i 



1 dxi dx2 

Xi X2 



5 1 



r 



Iz 



X\X2 



6{r] - ri{xipi,X2P2))Ff/h,{xi)Fj/^^^{x2). (29) 



^In fact, for the mass-distribution, the question can be fully answered up to two loops, given the recent 
work of pi. This is not so for the differential mass distribution, described in sec. 2.2, since the full two-loop 



results are not available at present. 



If we choose the kineinatical variable rj to be the rapidity y, 

r] = y = -ln{ , 30 

2 \qo-qLj 

we obtain for the partons 

y{xipi,X2P2) = 2 ^^ ( ^ ) ' (2^) 

so that a.t y = 0, 

''^'^ r,Q^y = 0)^a^lyJ:o(r,Q2), (32) 



rfV^ 



where _ 

WiL,{T,Q') = Y.CjA{a) J' dz[l + n{z,a)]-^^ (^^^^f^j , (33) 

with 

:Fyf{r/z) = Ff/,, (v^) F;/,, (v^) . (34) 

Similarly, choosing fj to be the Feynman parameter xp, 

f] = XF = 2qL/y/s, (35) 

where q^ is the longitudinal momentum of the virtual vector boson in the center-of-mass 
frame of the hadrons, we haveQ 

xf{xiPi, X2P2) =xx- X2 ■ (36) 

Therefore the resummed double-differential distribution at x^ = is 

^^^(r, Q\ x^ = 0) ^ 4(Q')W^i;=o(r, Q\ (37) 

where W^=o ^^ given by eq. (p3D, with the substitution 

J^=\rlz) -. r-=\Tlz) = ^=Fs/K. (^) F;/,, (v/^/l) . (38) 

2Jt/z ^ ^ ^ ^ 

Finite-order approximations can also be written down in a similar way: 



with 



a,b 



dz 



(40) 



•^In actual experiments, xi typically corresponds to the momentum fraction of a parton in the beam 
hadron while X2 corresponds to that of a parton in the target hadron. 

8 



where f] = y or xp and uj^^biz', a) is given by eqs. (|2^)-(|25|). 

We illustrate the accuracy of the approximation (|27|) in two ways. First, we compare in fig. 
1 the full 0{a) distribution (i.e., all 1-loop non-singular terms in the hard part are included) 
using eq. (|39| ) (with i = 1) and the exact distribution derived in |jl3|], in the central region. 
The kinematic parameters used are taken from the Fermilab E605 experiment |T^ and are 



described in section 4. The solid curve is the exact one-loop calculation while the dotted 
curve is the corresponding approximation. We see that the approximation is not perfect, 
with about a 10% difference between the two curves. This difference comes mostly from 
the quark-gluon channel. Hence, in all our subsequent calculations for double differential 
distributions, we used the exact results of |]13| for the quark-gluon channel alone. The 



remaining error is then of order 1%, which shows that the general approximation (eq. (|27|) ) 
for the channel containing the singular corrections is indeed very good. 

In most experiments, what is measured is not the fully integrated mass distribution 
da/dQ"^ but instead one that is partially integrated over either the rapidity y or xp in some 
kinematic range 7^(r). In order to compare with actual experimental data, we use the 
following approximation 

(S) =/ . ^^5:F^/ ,dvi7&- = ^B{Q')Wl^^^,ir,Q% (41) 

\dQ'^ J f](in(T) Jfi€n{T) dQ^drj Jfji^niT) dQ^drj ' ^ ' 

with d^a^ /dQ'^df] given by eq. (^7]) so that WY^^t^Ar^Q'^) is given by eq. (|33|) with the 
substitution 



^dxidx2 I 'r/z\ 

fi&l(T) 'Jo Xi X2 \ X1X2J 

x6{f] - r]{xipi,X2P2))Ff/h,{xi)Fj/,,^{x2). (42) 



^Jfir/z) - ^fj{r/z)^ l_^,^dni——6 



The second way, therefore, that we will test eq. (pTj) is by comparing eq. (^T]) in finite-order 
pQCD with an exact calculation at the end of this section. 

The integrals over xi, X2 in eq. (|^ are trivial to perform. Let us first consider the 
case of f/ = Xi? with 7^(t) the positive kinematic range (0, 1 — r), as is the case for most 
experiments. We will denote the corresponding function in eq. (|^ ) by jF^f ^ . Notice that 
upon performing the integrations over the parton fractions, the constraints Xi — X2 = xp and 
2^i^2 = "t/ z transform the positive kinematic range of xp into (0, 1 — t/ z), which makes the 
function JF^^^" depend on the scaling variable t/ z only. We obtain the result 



>o. / X /-^-^ dxp xp + ^xI + At/z \ [-xp+^xI + 4t/z 

^^^ ^^''^^L ^^TWz^'''\ ^ J ^"^" 

(43) 
On the other hand, ii t] = y and we want an analogous expression which holds for y > 0, the 
corresponding constraints are y = | ln(xi/x2) and X1X2 = t/z so that we get 

"^Note that xp > implies y > and vice-versa, and in fact using either xp or y for representing the 
differential mass distribution amounts to a change of integration variable, so that eqs. (|4^) and (|44| ) give 
identical results. We found it convenient to use eq. (|44| ) for calculating the y (or xp) > cross sections. 



J^^f(r/z) 



// 



imf 



dyFf/h^ 

Summarizing these results, we write 

'da"" 



-e^ ) F 



f/h2 



T 



dQ"^ J j)>o 



crimwL,{r,Q 



where 



Wlo^ 



r, Q2) = ^ C]A{a) / dz[l + 7^(2, a)] 



dz 



(44) 



(45) 



(46) 



Finite-order approximations take a similar form: 



\dQ^ 



fi>0 



^hmw;^o{r,Q 



(47) 



with Wy^Q^T^Q"^) given by eq. (0) but with J^^^\t/z) replaced by J^^^°{t/z), eq. (|g). 
In fig. 2 we test the approximation shown in eq. (41) by again comparing the resulting 

partially integrated distribution, ' '^'^ 



,„2 I , with the exact one using the results of [|T3l. The 

"'^ Jy>0 

calculations were done to one loop order, and the comparison is made within the context of 
the E772 experiment, with kinematics that we will describe in sec. 4. We see that in this 
case, the approximation works very well. 

In section 4 we will use these formulas in presenting our predictions and comparisons 
with various fixed-target experiments. 



3 The function A{a) 

We now turn to the coefficient function A{a) in eq. (|13D, which is closely connected to the 
Sudakov form factor. In [^, this function was given to 0{a) in the exponent as: 



Ai{a) = (1 + 2Cira)exp 



fc.(l.^-3 



(4J 



In this expression, the piece ^Cftt"^ of the exponent is identified with the absolute value 
squared of the ratio of the Sudakov form factor at timelike and spacelike momentum transfer, 

r(g2)/r(-Q2). 

To obtain an 0{a'^) form of A{a), we use eq. (3.8) of 0: 



In 



r(-g2 



IT 



a 



7^(1) + «2^(2) 



H ; — mv -\ ; — a 1 

262 262 V u 



(1) 1 (2) 

.2^_(^,lnzy - zy + 1) + 2|_(in^, _ ^, + 1) 
46| a ^ ^ Abl ^ ' 



(49) 



10 



where 



h -- 


12 ^ 


2 

12"^' 




V = 


= l+iTTb2a, 




li" = 


= 2Cp, 




7i? = 


= CaCf 


[?-H 


5 ^ 



(50) 



The 7]^- s are anomalous dimension coefficients. 

It is sufficient to consider only the real part of eq. 



Re 



In 



nQ' 



T{-Q') 



2bo 



Inr + 



TT 



G(2)sin^ 



-a 



(1) 1 (2) 

^^ (lnr-7rM«) + ^lnr, 



where 



1 + (7rb2a) 



The constant G^^^ is obtained from eq. (2.7) of 



461 a 

= arctan(7r62«) • 

i, 



461 



(51) 
(52) 

(53) 



while the constant G^'^^ can be found from their eqs. (4.1), (4.2), (4.3), (2.5) and (2.6)0: 



^(2) 



c'. 



1 _ !^(3) + ^^(2)' 
16 3^^ ^ 6 ^^ , 



+ CaCf 



432 12^^ ^ 4^^ ' 



-UfCf 



209 



:C(2) 



216 6 
Expanding eq. (|5ll) up to 0(q;^), one obtains the exponentiated result 



r(Q^ 



n-Q' 



a ^ 2 ^2 



exp -C^TT +-vr 



/233 2 
V~9~~ 3 



-TT^C^C^-^njC^ 



38 



(54) 



(55) 



The 0{a'^) exponentiated function A{a) can now be written as 



Aoia] 



1 + 2CFa + ha^ exp -Cf 






a 



-vr 
3 



+ 16" 



2 
/233 2 2\^ ^ 38 ^ 



(56) 



The constant 6 is determined by expanding eq. ( ^6]) in a and then comparing the result with 
the terms proportional to 5(1 — z) in the second-order result ^ \li eq. (^). This yields 



cl 



-23 



35 



15 



15 



720 



-vr 



7^^' + ifC(3) + ^ 



96 



+ CaCf 



215 175 



+ 



144 216 



-vr 



49 



17 
1440' 



-vr 



+nfCi 



rC(3) 



19 7 
72 ~ 54' 



rvr 



(57) 



^There is a misprint in eq. (4.3) of ||^. The right-hand side of the equation should read 



-4er 



(2) 
2 


'.^. 


" - 2.rf > 


'^\ 



(i^(i) + G(i))rw. 
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There is yet another way of treating the function A{a). The accuracy of the expansion of 
eq. (|5lD actually depends on the size of 7r62a- Given the value of 62 = 2.0833 at Uf = 4, one 
can easily check that 7rb2a is not really a very small number. Thus, instead of expanding in 
a, one can use the exact expression in eq. (|5lD to obtain the following expression for yl(a):[| 



A{a) = ( 1 + 2Ci.a + ba^j exp j -Cp 



3 



+ -— In r + 7tG^^> a^ 

ho r 



^(1) 1 ^(2) \ 



261 a ' ' 26, 

Note that at high energies, eq. (0) becomes a good enough approximation to eq. ( pT] ) so 
that the constant h in eq. (|58D is the same h that appears in eq. (|56|), given explicitly by 
eq. (13). 

In fig. 3 we have plotted the different forms Ai{a), ^2(0) and A{a) for nj = 4 and 
A = 0.25GeV. For consistency, we have used the two-loop form of a for ^2(0) and A{a) and 
a one- loop a for Ai{a). One sees that the following relation holds: A{a) > A2{a) > Ai{a) 
where ^4(0;) is given by eq. (p8|), A2{a) by eq. (|56D and Ai{a) by eq. (^8[) . As Q increases, all 



three functions become numerically close to 1, as expected from asymptotic freedom. On the 
other hand, at low values of Q all three functions have large numerical values, and diverge 
as Q approaches the confinement scale. In the mass-range of fixed-target experiments, all 
three forms of A are numerically very similar to one another. In the following section we 
will give predictions for the resummed dilepton cross section using only the exact, O(a^) 
function A{a). 

4 Comparisons with experiment 

4.1 Generalities 

In this section, we present our numerical calculations and compare them with results from 
several fixed-target experiments. To be consistent with the two-loop running coupling used 
in the calculation of the resummed exponent , we use the following two-loop form of a in 
terms of A: 

(n^^ = "^(^') = 1 _ 63 ln(ln(QVA2)) 

^^^" TT 62ln(QVA2) 63 ln2(Q2/A2) ' ^^^^ 



with 



62 = {UCa - 2nf)/12, 63 = (34Cl - {IOCa + 6Cp)r2/)/48. (60) 

For parton distributions, we used the recent CTEQ2D |jl4] and MRSD- |1^ parton sets 



that are based on global fits. In the hard part, we used the A„ values provided with the 
appropriate parton distribution. The parameter nj is either 4 or 5 depending on the value 
of Q at which the parton density is evaluated. 

We note here that in all the discussions to follow, by resummed cross section, we mean 
the resummed cross section obtained from eq. (p2D or (|45|) plus the contributions of all the 



^ In the following, it is this "exact" expression that wc will denote by A{a) and use in the numerical 
computations. 
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non-singular terms in the 1-loop hard part, ujI}, {z^ a). Specifically, these are ctUjg (eq. (|2^)) 
and a^{—Az — 6 + 2(1 — z) ln(l — z) — 41n(l — z)), derivable from aUrj- of eq. 

4.2 The NA3 (1980) Experiment 



This experiment ||T^ involved the hadronic production of dimuons from the interaction of 
a proton- {^/s = 19.4GeV) and anti-proton- beam {^/s = 16.8GeV) with a platinum (Pt) 

target. In fig. 4a we have plotted the data points together with the one-loop, two-loop 

3 da 
and resummed forms of the mass distribution M,, 



using the CTEQ2D (A4 

xp>0 



0.235GeV) distribution. For both p— and p— on Pt, all three curves exhibit acceptable 
agreement with data. 

We show in fig. 4b the resulting curves using the MRSD— (A4 = 0.23GeV) parton set. 
For p on Pt, the curves are almost identical to the corresponding ones in fig. 4a, while the 
p on Pt curves are slightly higher than their CTEQ2D counterparts. These suggest that 
CTEQ2D and MRSD— have almost the same valence quarks but slightly different sea quark 
parametrizations. 

It should be noted here that the size of the resummed cross sections is largely determined 
by the function A{a). In fig. 5, by looking at the dotted (Born) and dashed (resummed but 
with A{a) replaced by 1) curves, we note that the contribution of the plus distributions 
steadily increases with i/r, from only a few percent in the range ^/T < 0.3, up to about 26% 
at y/r c^ 0.5. Comparing the dashed curve with the solid one, we can also determine that 
the plus distributions contribute from only a few percent to about 11%, in the same regions 
of \/r, to the resummed cross section. One can also check that, in the mass range of this 
plot, at least 80% of the resummed cross section comes from A{a) x Born cross section. 

4.3 The NA3 (1985) Experiment 

In this experiment [jl^, dimuon events, at ^/s = 21 A GeV, produced in the collisions of 
protons with a platinum (Pt) target, and within the range xp > ^ were analyzed. 

In fig. 6 we present our resummed prediction for the mass distribution Mf^ —— — , 

^^dM^^ xp>o 

using the CTEQ2D (A4 = 0.235GeV) parton set. The resummed curve is no more than 23% 

above the data, and the 1-loop prediction fits the data best. 

4.4 The E537 Experiment 

This Fermilab experiment [|l^ involved collisions of an antiproton beam with a tungsten 
(W) target at ^/s = 15.4 GeV. Again, dimuon production events which fell in the positive 
xp-YSiJige were selected. 

In fig. 7, we have plotted the data points and the mass distribution M^^ 



'dM^^ 



usmg 

xp>0 



the CTEQ2D parton distribution. Here, the resummed curve has the best fit. The 2-loop 
curve also has a good fit but with a shght underestimate. 
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We have also plotted the differential mass distribution M^ 



,,., using the 1-, 2- 

dMdxF xp=o 
loop and resummed hard parts in fig. 8. The resummed curve has an excellent fit, while 

both finite order curves underestimate the data. Again, as in the NA3(1980) p on Pt case, 

the error bars are so large that all three curves are within them. 

4.5 The E605 Experiment 

In this Fermilab experiment |jl9|], a proton beam was incident on a copper (Cu) target 
at -yi = 38.8 GeV. We show in figs. 9a and 9b the data points with the 1-loop, 2-loop 

using the CTEQ2D and MRSD- distributions, 

y=0 



and resummed predictions for s- 



dy^dy 
respectively. 

E605 data together with a 1-loop hard part as a theoretical input, were used for the 
CTEQ global fit, so, as expected, the 1-loop curve in fig. 9 gives an excellent fit to the data 
over the whole range of r. The resummed cross section is about 18 to 14% bigger than the 
data, but within experimental errors, the difference decreasing on average as r increases. 

In fig. 9b, the 1-loop curve also has the best fit to the data but with some overestimate. 
The 2-loop and resummed curves overestimate the data, and are about 10% bigger than 
their CTEQ2D counterparts in fig. 9a. 

4.6 The E772 Experiment 



This recent Fermilab experiment |^ focused on interactions of a proton beam with various 
targets, to obtain, among other things, the nuclear dependence of the resulting spectrum. 
Here we will look at p-^H cross sections at ^/s = 38.8 GeV and xp > 0- 



In fig. 10, we plot the mass distribution Mf 



""dM^^ 



using the CTEQ 2D parton dis- 

xp>0 



tributions. In the low mass region (.1 < ^/T < .24), the resummed curve has the best fit 
to the data. In the high mass region (.1 < y/r < .24), the 1-loop curve has the best fit. 
However, the resummed and 2-loop curves are no more than 16% bigger than the 1-loop 
curve, and within experimental errors. 

We show in fig. 11 the distribution Mf, —-— — ; — using CTEQ2D parton distribu- 

^^dM^^dxp xp=o 
tions. In the high mass region, all three curves pass through almost all of the error bars. The 

resummed and 2-loop curves, which in this case are almost identical, have the best fit. In the 

low mass region, the 2-loop cross section gives a good fit, with some slight underestimate. 

The resummed curve overshoots the data by 18% or less. 

4.7 A'- factors 

We have also calculated a set of K-factors for the kinematics of various experiments using 
eqs. ([I2D, (|17D , (|2l| ) and ([19|) . Note that the denominator is simply the Leading Logarithm 
Approximation (LLA) cross section. For figs. 12 and 14, the K factors were calculated over 
the full rapidity range, while for fig. 13, the K^s were only calculated over positive rapidity. 
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These quantities express the size of the resummed radiative corrections relative to lowest 
order. They turn out to be relatively insensitive to the parton sets used, and allow us to 
compare the importance of corrections for the various experiments. 

In fig. 12, we plot Ki, K2 and the resummed K factor, using parameters from the E537 
experiment described above and the CTEQ2D parton distributions. Ki is the z-th loop K- 
factor obtained using eq. (pll). To better illustrate the different contributions to K, we show 
two other curves in the same plot: the dot-dashed and dotted curves are resummed curves 
but with A{a) (eq. (|58|)) replaced by Ai{a) (eq. (|i8D) and 1, respectively. 

We see that, in the mass range covered by the plot, most of the resummed K comes 
from the A{a) term, as the plus distributions in this case account from only a few percent to 
about 26% of the total K. This agrees with an earlier assessment given in the cross section 
discussions. Comparing the solid and dot-dashed curves, one also observes that the next- 
to- leading terms in the exponent of A{a) yield from 3 to 7% contribution to the resummed 
cross sections. 

In fig. 13, we present our theoretical resummed X-factor, along with Ki and K2 for the 
NA3 (1985) experiment, using two different parton sets: CTEQ2D and MRSD— . We notice 
that the resulting curves are in very good agreement with one another. 

In fig. 14, we present the i^-factors, using the CTEQ2D parton distributions for the 
E772 experiment. What is important to note here is that as one goes to large mass values, 
the difference between the resummed and the 2-loop i^-factor decreases to about a few 
percent. Moreover, whereas Ki and K2 are noticeably increasing in this region and beyond, 
K is relatively flat or very smoothly decreasing. This highlights the "well-behaved" nature 
of principal value resummation, for it is precisely in this region where the logarithms in 
perturbative calculations are large. 

4.8 Theoretical Uncertainty 

In this section, we quantify the theoretical uncertainty in our resummed calculations, as 
well as the goodness-of-fit of both resummed and finite-order cross sections. For purposes of 
discussion, we focus on data from the E605 and E537 experiments. 

We compare in fig. 15 the resummed cross sections calculated using CTEQ2D and 
MRSD— parton distributions with data from E605. We see that the 2 curves are within 
10% of each other. There is also an estimated uncertainty of about 7% due to contributions 
of non-singular terms that were not included in the resummed hard part (see discussion in 
section 5). Thus, we consider a good estimate of the theoretical uncertainty in the resummed 
predictions to be 17%. 

In table 1, we present the x^ per degree of freedom of the calculated cross sections relative 
to data from E605. As expected, the 1-loop cross sections calculated with global parton sets 
have the best x^ values. However, one also notes that the x^'s of the resummed cross sections 
are not that large compared with the best values. [] 

^In M, we discuss the possibility that this difference is related to nonperturbative effects. 
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TABLE 1 



X^ per degree of freedom (E605) 



Parton set Resummed 1-loop 2-loop 



CTEQ2D 



5.0 



1.3 



2.1 



MRSD- 



11.6 



2.6 



5.9 



Similarly, we present in table 2 the analogous x^ values for the experiment E537, using 
the data points shown in our fig. 8. In this case, the resummed cross sections have the best 
X^ values, while the finite-order cross sections have x^'s not more than 1. This supports our 
earlier claim that for antiprotons on fixed targets the resummed cross section has the best 
fit. 

TABLE 2 



X per degree of freedom (E537) 



Parton set Resummed 1-loop 2-loop 



CTEQ2D 



0.05 



0.38 0.84 



MRSD- 



0.14 



0.46 0.93 



5 Discussion and Conclusions 

In this paper, we have used a previously calculated hard part [@] for dilepton production, 
which resums large threshold corrections to all orders in pQCD via the method of principal 
value resummation 0, to calculate resummed DY distributions. We have also derived the 
coefficient function A{a), which involves the resummation of large, 2;-independent Sudakov 
terms, to next-to-leading order. We then made comparisons of several distributions with 
data from fixed-target experiments. In as much as the available data cover only the positive 
rapidity region, we used the approximation of |]12[ in our calculations. All calculations have 
been made in the DIS scheme. To be entirely consistent, one should use parton distributions 



resulting purely from fits to Deep Inelastic Scattering data. Existing DIS fits (like DFLM |^ 
and EHLQ |23]), however, had been constructed from relatively old data. On the other hand 



the global fits (like CTEQ2D and MRSD— ) that we have used, typically use only 1-loop hard 
parts for input in hadron-hadron scattering, for example when Drell-Yan data is included, 
but include the most recent DIS data. If higher order corrections are important, their effects 
will be included in the globally fitted distributions, and using them in a resummed calculation 
would involve some double-counting and thus lead to overestimates. 

We now summarize our results in the light of these issues. For all fixed-target experiments 
we have examined, the resummed cross sections exhibit good to excellent agreement with 
the data. For antiprotons on fixed targets, we find in general that the resummed predictions 
agree better than the finite-order ones. For proton beams the finite-order predictions agree 
better in most cases, but in those cases the resummed ones are above the data by no more 
than 10 — 23% and even then, all predictions are within experimental errors. From the above, 
we can conclude that principal value resummation has successfully "tamed" the divergences 
in perturbation theory associated with the resummation of plus-distributions. We find that 
even if large corrections exist at the 1- and 2-loop order, principal value resummation yields 
cross sections that are finite and well-behaved. In addition, the quantitative structure of 
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the large threshold corrections in principal value resummation correctly accounts for the 
observed data. 

It is important here to summarize what these results represent concerning our knowledge 
of perturbative QCD. Our resummed cross sections were calculated including all the radiation 
phase space down to threshold. This includes the perturbative regime, which provides the 
bulk of the effect, and the non-perturbative regime which we treat via the principal value 
resummation model. 

The fact that, using global parton sets, the resummed cross section slightly overshoots 
most proton beam data, needs to be disentangled: 

It would be useful if a parton set fitted only to the most recent DIS data becomes available, 
so that one could repeat our resummed calculations with this set. If, however, the antiquark 
fits are so "robust" that the resulting cross sections do not change significantly then the 
following issues, which could possibly explain the overestimate, need to be confronted. First, 
the cross section could be sensitive to higher twist. The physical higher-twist would be 
slightly different than in principal value resummation, and should be parametrized more 
generally. This issue is studied in more detail in [^ . 



Secondly, the factor A{a) accounts for a large part of the resummed cross section. It is 
important to note that only the exponentiation of Sudakov terms in A{a) is fully understood 



at this moment. Since the exponent of eq. (|58|) contains some finite non-Sudakov terms as 



well, following pj] , a better understanding of the resummation of these terms in this function 
is also needed to eliminate some of the theoretical uncertainty. It might be that inclusion of 
such terms might make A{a) smaller, so that the resummed cross section is also reduced. 

Another issue is the size of the non-singular terms in the hard part at higher orders. For 
the case of the 2-loop order, the size of the contributions of these non-singular terms has 



been studied. From |TT|, we note that in the kinematic region we have considered, the non- 



singular terms in the 2-loop hard part yield from 5 to 10% negative contributions to the cross 
section. In the resummed case, those terms will lead up to a 7% negative contribution.^ If 
such "finite" terms continue to be sizeable at still higher orders, however, they might further 
reduce the effects of resummation, which might explain why 1- or 2-loop cross sections that 
use global parton sets can describe the data quite well. 

All of these issues are correlated in the differences of our predictions using different parton 
sets, and the 2-loop reduction of the non-singular terms. Our uncertainty estimate of 17% 
is in fact a conservative measure of these factors and clearing-up these issues corresponds to 
narrowing down this theoretical uncertainty. 

In future work we hope to extend these calculations in the context of collider energies, 
both for the Drell-Yan continuum and for vector boson production. Among the most inter- 
esting other reactions sharing many common features with dilepton production, and where 
principal value resummation has been advantageously applied, is top-quark production ||25|| . 
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APPENDIX A Some Notes on the Numerical Calculations 

The explicit expression used to calculate the resummed cross section in this paper is 

da sr^ 2 f^ 7 A/ \t -Ml \T ^ f^ffi^/^) 



dQ^ 



^bY.^)1 dzA{a)[l + n{z,a)]^{-^f^^\ , (61) 



where as is given by eq. ([Zp, A{a) by eq. (^BD, 7Y(z, a) by eq. ( ]T^ ) and the parton luminosity 
by eq. (^. The factor involving the luminosity can be expressed as an integral over the 
rapidity y g: 

d (J^fj{T/z)\ _ I fy^ , f^ d d 



dz\ z J 1z 



y2 



l[ dy(2 + x^— + X2—yf/HA^uQ')Ff/,^ix2,Q'), (62) 



where Xi = a/^g^, X2 = J^e ^ and y^ = ±ilog(^). 
Note that eq. (^) may now be re-written as 

da dao v^ 2 /"^ r -i// \f ^ \ P^ , fc^ ^ d 



AW^ + ^(aK E 4 />--«(--.") (-^) /; ''K^ + ..A + ,,. 



^2 



yiFf/hAxi.Q)Ff/h^{x2,Q), (63) 

with j^ as in eq. (p|). With the integral definition of 'H{z,a), eq. (|18D, one has to per- 
form at least four integrations to compute da/dQ^ for one value of Q^ Straightforward 
implementation of eq. (pT| ) is thus very time-consuming from a numerical point of view. 

In order to shorten significantly the processing time requirements, an intermediate step 
was performed: the quantity Ti.{z, a) was first interpolated over a range of z and a values. 
The interpolation is essentially a two parameter {z, a) fit to a set of H{z, a) values calculated 
over the box contour for particular values of z and a. But because E (and thus Pi) Q depends 
on parameters which in turn depend on n/, there is an implicit dependence of 7i{z,a) on 
Uf. Furthermore, note that one can have a leading or next-to-leading T-C{z, a) depending on 
whether E = E^ or E = El + E^l is used. For all the fixed target cross section calculations, 
ri/ = 4 and E = El + Ej^l were used. 

Armed with the interpolating function, the calculation of eq. ( |63D becomes a little sim- 
pler (for now, only a 2-fold integral is involved) and more importantly, a lot faster from a 
numerical point of view. 

The interpolating function, which we will now denote by 'Hi{z, a), is valid for 
.01 < z < 1 - 10"^ and 4.756 x 10"^ < a < 7.197 x 10"^ The choice of the boundaries for 
z is essentially machine-dependent: outside this range in z, the program ran into round-off 
errors. Hence, when Tii^z^a) is used in eq. (^), the upper z-limit is taken to be 1 — 10^^. 
As mentioned in our previous paper 0, for fixed a, the exact Tii^z^a) is bounded, reaching 
a certain limiting value before 2 — ;> 1 — 10~^. The contribution from the region z > 1 — 10~^ 
is thus expected to be very small and may be neglected. 

^ Aside from the z, z' and y integrals, further integrations are required to calculate E and Pi. 
^°In calculating the set of 7i(z, a) values, E was calculated over a box contour g while Pi was calculated 
numerically: Pi{z) = (1 — z) ^^ '^ ^^~ ' for some small number h. 

-'^-'^The interpolations were obtained using the B spline routines in the Fortran 77 IMSL package. 
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For the finite order calculations, instead of eq. (]6T|), (|62D, we use 

^Fa/hA^l^Q)Fb/h2{x2,Q), (64) 

with C^ given in section 2.1 and ul^jj{z',a) given by eqs. (|2^)-(p5D. No interpolations were 
necessary for the finite-order hard-parts. 

For the exact calculation of 
expression for the qq channel is 



da 



(1) 



dQ'MyW'- 



=0, we used the results given in |T^. The explicit 



da, 



(1) 



dQ'^dy 



_ 167ra|a^ 2 v"^j 

=0" 27 g2, 1.^/1.//, 



//' 



(65) 



with 






1 + Ivr^ - I In TT-^ — T + 2 In jr^ In j-^ 

3 2 (l-xi)(l-x2) (l-^'i) (l-a;2 



P./(xi,X2,(5), 



^ff ~ 2 Jxi "''lU^(ti-a;i)+ "^ X2(ti+xi) "^ 2 (ti-xi 



--^-^)^//('^i'^2,Q), 



^// - 2 /i ^^2 (^i^(4^!^2)+ In ^^(i2+^2) + I (t2-l2)+ " I ~ ^ J%(^1' ^2, Q), 

^// = /ii ^^1 /i ^^2 [(i^_^^)(*j/42)]^ ^//(^l' ^2, <5), 
% = Li ^^1 L2 dt2t^t^(^tJ^+ll^y2Pff{tl,t2, Q), 

where xi = X2 = ^r, Pff{ti,t2,Q) = Ff/h^{ti,Q^)Ff/h.,it2,Q^) and 

^'^(ti t2'- (^ + '^i^2)(r2 + (tit2)') 



(66) 



(67) 



(tlt2)2(tl+Xi)(t2+X2)' 

The sum in eq. (^5|) extends over all active quark and antiquark flavors. Eqs. (|65|) and (|66|) 
were obtained using eqs. (2.1), (2.3), (2.6)-(2.8) of [0. The -|- distributions are evaluated 
using the prescription given by their eqs. (2.9) and (2.10). 

Similarly, for the quark-gluon channel, we used eqs. (2.13)-(2.17) of |T^ to obtain 



HJ 



dQ'^dy 



y=o 



27r aja 



f i=i 



(6^ 



with 



rrl _ fl ^f hrj + {ti-xif , 2x1(1-3:2) I 1 3xi(ti-xi) \ pC (f „ -Q\ 



Ixi dtl /^2 dt2-lj^z^^Pgf{ti,t2, Q), 



Ulf = I^,dhJ^^dt2H^{h,t2)P^f{tl,t2,Q), 






gf^ 

xl + {t2-X2f 1^ 2X2(1-Xl) 






'"^^*^-"^Mif;(t2,Xi,Q), 



^1/ = /4^t2/irft,fi|^pC'(,^^,^^g)^ 
^,'/ = I^,dt2j^^dhH^it2,h)Pg{t2,h,Q), 
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(69) 



where 






X2{j+tit2){T'^ + (r-tlt2f) 
i-{i^{tlX2+t2Xi)(t2+X2) ' 

Lt2)(*li2^1+T(tiX2H 
{tit2)'^{tiX2+t2Xi)^ 



PSit2,h,Q) = Fg/H2it2,Q')Ff/h,ih,Q^) 



For the exact calculation of 



da^ 



dQ^ 



y>o 



, we used 



da, 



(1) 



dQ^ 



y>o 27 Q'^s 






a _ 



with 



(70) 



(71) 



ci_ 

C2_ 



1 fl_ dxi 
2 



1 rl dxi rl jf ( tj+xl , 2xx{xi-t) . 3 1 2 _ 3x1 \ p _/. / f)\ 
2JV^ — Jxi "tl ytl(h~x,)+ ^^ r(ti+a;i) + 2(ti-a;i)+ ~ IT IfJ^ff^^^^ ^/ ^1' ^''' 

--t)%(^/^2,t2,Q), 



^2~^^2 



jj^ 2x2(x2-t) _|_ 3 



rVT 0X2 n A-j- 
2Jt X2 ■'X2'^^^\ll{t2-X2)+''''^ T{t2+X2) ' 2(t2-a;2)+ t2 



^ff - /v^ ir -f^i ^^^ -fr/xi dt2 yi^t^-x^)(t2-l/x^)]+Pff^^^^ ^2, <5), 
qS _ 



/v^ ^ /xi d'^l /r/xi ^^2t^t2(tirAri+i2Xi)^^//('^l' ^2, Q)- 



(72) 



Eqs. ( |7lD and ( [72D were obtained using eqs. (2.6) and (C.l) of [|T3l. Note that the integration 
limits for Xi, X2 follow from the requirement that |/ > or equivalently, xp = Xi — X2> 0. 

In a similar way, starting from their eqs. (2.17) and (C.l), the following formula is ob- 
tained for the quark-gluon channel: 






dQ^ 



y>o i) (^ s J. j^-^ 



(73) 



with 



Kf 



rl dxi rl J. ( xl+(ti-x^Y 1^ 2xi(a'i-r) , 1 3a'i(ti-xi) \ pC /y. ^/™ r)^ 
Jv^ ~ Jxi "'^l (^ 2tf ^^ r(ti+xi) + 2*7 i? J^gA^^' ^/^l' ^''' 



^// = //'^/.V2rf^i/i2^^2SS^p,^,(ti,t2,g) 



9/ 



f!/T ^ £, *i /,V, dhH'={tuh){f,(h,h, Q), 



T/4 _ rV? dx2 fl Jj- / x^+(t2-X2)'^ 
^9f ~ Jr — Jx2"^2| 



2K 



1,^ 2x2(a:2-r) , J_ 

r(t2+a;2) "^ 2*2 

K'f = Sy-ltU.dt2tdti^^^PSit2M.Ql 



'-^^^^^^^)P^;{h,T/x2,Q), 



Iv-r ^ Ir/x, dh t dtlH^{t2, tl)P^;{t2, h, Q) . (74) 

All of the various parton distributions used were taken from the CERN PDFLIB package. 



20 



FIGURE CAPTIONS 



Figure 1. Exact and approximate s- 



rf2 



a 



to 0{a). 



dy/rdy y=o 

Solid=exact calculation using results of Kubar et al in ITsI ; 
Dotted=calculated using Laenen-Sterman approximation. 



o da 
Figure 2. Exact and approximate M,,„-— -— 



MM 



xp>0 



to 0{a) 



Solid=exact calculation using results of Kubar et al in |T3| 
Dotted=calculated using Laenen-Sterman approximation. 

Figure 3. Different forms of A{a). 

So\id= A{a); Dashed=A2(a); Dotted=Ai(a). 



Figure 4. NA3(1980): M; 



da 



MM 



dM, 



/ifi 



XF>0 



(a) CTEQ2D. Solid=resummed; Long Daslied=2-loop; Short Dashed=l-loop. 

(b) Same as (a) but for MRS D-. 

da 



Figure 5. Contributions to M 



fifj. 



dM, 



tlfM 



XF>0 



Solid=total resummed cross section; Dashed=resummed but with A{a) replaced by 1; 
Dotted=Born cross section. 



Figure 6. NA3(1985): M^ 



da 



IMfl 



dM, 



/IfM 



Xp>0 



using CTEQ2D partons. 



Solid=resummed; Long Dashed=2-loop; Short Dashed=l-loop. 
da 



Figure 7. E537: M. 



iifi 



dM, 



Hfj. 



xp>0 



using CTEQ2D partons. 



Solid=resummed; Long Dashed=2-loop; Short Dashed=l-loop. 



Figure 8. E537: M-" 



dMdxi 



xp=0 



using CTEQ2D partons. 



Solid=resummed; Long Dashed=2-loop; Short Dashed=l-loop. 



Figure 9. E605: s 



d^/Tdy 



y=0 



Solid=resummed; Long Dashed=2-loop; Short Dashed=l-loop. 

(a) CTEQ2D. 

(b) MRSD-. 

da 



Figure 10. E772: M. 



fi/j. 



dM, 



fifi 



xp>0 



using CTEQ2D partons. 



Solid=resummed; Long Dashed=2-loop; Short Dashed=l-loop. 
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Figure 11. E772: M^ ^ 



using CTEQ2D partons. 

xp=0 



dMdxp 
Solid=resuninied; Long Dashed=2-loop; Short Dashed=l-loop. 

Figure 12. K factors for E537 p-W using CTEQ2D. 

Solid=resumnied; Long Dashed=2-loop; Short Dashed=l-loop; 
Dot-Dashed=resummed with A{a) replaced by Ai{a)\ 
Dotted=resunimed with A{a) replaced by 1. 

Figure 13. NA3(1985) theoretical i^-factors for different parton sets. 
Dashed (descending order): Same as above using MRSD-; 
Dotted (descending order): Same as above using CTEQ2D. 

Figure 14. K-factors for E772 p-^H using CTEQ2D. 

Solid=resummed; Long Dashed=2-loop; Short Dashed=l-loop. 

Figure 15. Resummed curves for E605 using global fits. 

Upper solid curve=MRSD-; Lower solid curve=CTEQ2D. 

References 

[1] G. Altarelh, R.K. Ellis and G. Martinelh, Nucl. Phys. B157 (1979) 461; 
B. Humpert and W.L. van Neerven, Phys. Lett. B84 (1979) 327; 
J. Kubar-Andre and F.E. Paige, Phys. Rev. D19 (1979) 221; 

K. Harada, T. Kaneko and N. Sakai, Nucl. Phys. B155 (1979) 169; B165 (1980) 545 (E); 
T. Matsuura and W. L. van Neerven, Z. Phys. C38 (1988) 623; 
T. Matsuura, S.C. van der Marck and W.L. van Neerven, Phys. Lett. B211 (1988) 171. 

[2] L. Alvero and H. Contopanagos, Nucl. Phys. B436 (1995) 184. 

[3] H. Contopanagos and G. Sterman, Nucl. Phys. B419 (1994) 77. 

[4] D. Appell, P. Mackenzie and G. Sterman, Nucl. Phys. B309 (1988) 259. 

[5] T. Matsuura, S.C. van der Marck and W.L. van Neerven, Nucl. Phys. B319 (1989) 570. 

[6] L. Magnea and G. Sterman, Phys. Rev. D42 (1990) 4222. 

[7] R. Hamberg, W.L. van Neerven and T. Matsuura, Nucl. Phys. B359 (1991) 343. 

[8] CTEQ collaboration (G. Sterman, editor) Handbook of perturbative QCD, version 1.0 
(April 1993). 

[9] H. Contopanagos and G. Sterman, Nucl. Phys. B400 (1993) 211. 

[10] W.L. van Neerven and E.B. Zijlstra, Nucl. Phys. B382 (1992) 11. 

[11] P.J. Rijken and W.L. van Neerven, Phys. Rev. D51 (1995) 44. 



22 



[12] E. Laenen and G. Sterman, The Fermilab Meeting, DPF'92, Vol. 2, 987 (C. Albright et 
al., editors). 

[13] J. Kubar, M. Le Bellac, J.L. Meunier and G. Plant, Nncl. Phys. B175 (1980) 251. 

[14] CTEQ Collaboration, unpublished; H.L. Lai et al., MSU-HEP-41024/CTEQ-404 (Oc- 
tober 1994). 

[15] A.D. Martin, W.J. Stirling and R.G. Roberts, Phys. Rev. D47 (1993) 867; Phys. Lett. 
306B (1993) 145. 

[16] J. Badier et al., Phys. Lett. B96 (1980) 422. 

[17] J. Badier et al, Z. Phys. C26 (1985) 489. 

[18] E. Anassontzis et al, Phys. Rev. D38 (1988) 1377. 

[19] G. Moreno et al, Phys. Rev. D43 (1991) 2815. 

[20] P.L. McGaughey et al. (E772 Collaboration), Phys. Rev. D50 (1994) 3038. 

[21] L. Alvero, ITP-SB-94-67 (December 1994), |hep-ph 9412335| . 

[22] M. Diemoz, F. Ferroni, E. Longo and G. Martinelh, Z. Phys. C39 (1988) 21. 

[23] E. Eichten, I. Hinchliffe, K. Lane and C. Quigg, Rev. Mod. Phys. 56 (1984) 579. 

[24] G. Sterman, Nucl. Phys. B281 (1987) 310; 

S. Catani and L. Trentadue, Nucl. Phys. B327 (1989) 323; B353 (1991) 183. 

[25] E. Berger and H. Contopanagos, ANL-HEP-PR-95-31 (July 1995). 



23 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



> 

O 

c 



>> 

>^ 

"O 



CM 



T3 



CM 



"D 
CO 



p Cu (38.8 GeV) 



CTEQ2D (A4=0.235GeV) 



100 




.18 



.22 



^V2 

Fig. 1 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig3-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig4-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



<N 



> 

o 

C 



o 

A 



I 



b 

T3 



I 



ro 



10 



.1 



p ^H (38.8 GeV) 



CTEQ2D (A4=.235GeV) 




.15 



T 



1/2 



.2 



Fig. 2 



This figure "fig3-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig4-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig3-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



(A=0.25 GeV, n^=A) 



1 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ r 











J \ \ \ \ \ \ \ \ \ \ L 







20 



40 



60 



80 



100 



M (GeV) 
Fig. 3 



This figure "fig4-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig2-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig3-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig4-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



<N 



> 

o 



o 

A 



I 



b 

T3 



I 



ro 



P,pPt 



10 



1 \ r 



1 \ r 



CTEQ2D (A4=.235GeV) 



pPt(16.8GeV) 



1 \ r 



NA3 




.2 



.3 



.4 

1/2 
T 

Fig. 4a 



.5 



.6 



<N 



> 

o 



o 

A 



I 



b 

T3 



I 



ro 



P,pPt 



10 



1 \ r 



1 \ r 



MRSD- (A4=.23GeV) 



pPt(16.8GeV) 



1 \ r 



NA3 




J \ L 



.2 



.3 



.4 

1/2 
T 

Fig. 4b 



.5 



.6 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig2-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig3-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig4-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



<N 



> 

o 



o 

A 



I 



b 

T3 



I 



ro 



10 



.2 



p,pPt 



CTEQ2D (A4=.235GeV) 



1 


III 1 1 1 1 1 1 1 1 


1 1 1 


1 


- 




NA3 


— 


- 


pPt(16.8GeV) 




- 


- 


-~> ^-v^ 




— 


— 


'■■■- ^\ "•~"~ ^\> 




— 


- 


''-,^^ \\ "■■■■>-. 




— 




"''--. \x ''■■^" 








■>. \. pPt(19.4GeV) 

\ \ 






1 


III 1 1 1 1 '''-."^M 1 1 1 


1 1 1 


1 



.3 



.4 

1/2 
T 

Fig. 5 



.5 



.6 



This figure "figl-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig2-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig3-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



This figure "fig4-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411294v2 



<N 



> 

o 

C 



o 

A 



I 



b 

T3 



I 



ro 



10 



p Pt (27.4 GeV) 



CTEQ2D (A4=.235GeV) 



- 


- 


^^ NA3 


— 


^^"-^' - ^ 


^ 






L^:f 1 


^ 




^J 


L.^.K -. 


^^ 




J 


L-.:i^n 


-^^ 




"■ 


s, ^^ 






^-f ^T 






T" J 


^ 


^ 






^^* ^ 


^^ 






^ ^ f \ -1^^ 




- 


- 


^ 


^ 


^ 




_ 


L"-H, ^ 


^^^ 






"^ -V. •" ^-~\ 






^ 


- ^^ 


^^ 




^ J 


L^^-,k\-^ 


^\^ 






^f: 4^ 




-L 


-'C rs 








L 


^ 


\ 


^^v^ 






"^-:k ^ 


^ 






L 




^ \ 


^\^ 











^ 


^ 




- — ,— 





_ 


_ 


"«•-. -1 










_ -^ 


^ O \ 


\^ 




— 






^^-^ CK 










^ ^ 


^ 1 


- 




J 






^ \ 


^\^ 








\ 


~\^ 


\ 


— 


_ 


(» " ^ ^ 


\ 






( 1 




- 


- 





.2 



T 



.1/2 



.3 



Fig. 6 



10 



<N 



> 

o 

C 



o 

A 



I 



b 

T3 



I 



ro 



pW(15.4GeV) 



CTEQ2D (A4=.235GeV) 





1 


1 1 1 


II 1 1 1 1 

E537 


_ 


( 


;^^ 


: 


- 




" " ^ <k ^^^ -r- — 


— 




^ 


^ ^ ^ \ 


^\^ — 1 — 








^ ^ \ 






— 






\ 














"-^ 




— 










"- \ \ ~ 


— 










-^^ \ x. 


- 












— 














1 


1 1 1 






1 1 1 



.3 



^V2 

Fig. 7 



.4 



.5 



pW(15.4GeV) 



CTEQ2D (A4=0.235GeV) 



> 

0) 

O 

C 

o 

II 
li. 

X 



X 

"D 

:^ 

b 



-^ 10 



( 


1 
1 


1 1 
1 1 


1 
1 


1 

E537 

» _ 
1 



^V2 

Fig. 8 



.4 



> 

O 

c 



>> 

>^ 

"O 



100 



(N 



T3 



CM 



"D 
CO 



10 



p Cu (38.8 GeV) 



CTEQ2D (A4=0.255GeV) 




.2 



r 



1/2 



.3 



Fig. 9a 



> 

O 

c 



>> 

>^ 

"O 



CM 



T3 



CM 



"D 
CO 



100 



10 



p Cu (38.8 GeV) 



MRSD- (A4=0.25GeV) 




.2 



r 



1/2 



.3 



Fig. 9b 



<N 



> 

o 

C 



I 



b 

T3 



I 



ro 



10 



^^ 1 



.1 



p ^H (38.8 GeV) 



CTEQ2D (A4=.235GeV) 




.2 



^V2 

Fig. 10 



.3 



.4 



> 

CD 

o 

C 

o 

II 
li. 

X 

lT 
X 
"O 

"O 

b 



n 



10 



.1 



p ^H (38.8 GeV) 



CTEQ2D (A4=0.255GeV) 



- 


1 1 


II II 


II III 


- 




"fel-^ 




E772 




— 


- 




— 




~tc^ 


k 










^^^^\ 










^'••^>>\ 






— 






— 


I 




1"%-^ 




I 


- 




1^%. 




— 


- 




~r^ 




— 


— 






( ) 


— 








- %^, 


— 










^ C^\, 


t 




— 






^ t5^ 


r- 


— 


- 








■J 


K^°V^ 


— 


— 


- 








X 
■^ 




t 


— 


— 








( 


1 


^^^^ 


N4. 


- 


1 1 


II II 


1 1 1 


1 


1 


^- 



.2 



T 



.1/2 



.3 



Fig. 11 



pW(15.4GeV) 



(CTEQ2D) 



1 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ r 

E537 







J \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ L 



.3 



.4 



T 



.1/2 



.5 



.6 



Fig. 12 



p Pt (27.4 GeV) 







.2 



T 



.1/2 



(n,=4) 




.3 



Fig. 13 







.1 



p ^H (38.8 GeV) 



(CTEQ2D) 




.2 



T 



.1/2 



.3 



.4 



Fig. 14 



p Cu (38.8 GeV) 



> 

^ 100 

c 



>> 

>^ 

"O 



CM 



T3 



CM 



"D 
CO 



10 




.2 



r 



1/2 



.3 



Fig. 15 



